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Abitrad 

The use of Linear Quadratic Regulator (LQR) control synthe- 
sis techniques implies the availability of full state feedback. 
For vibration control of structures, usually only a limited 
number of states are measured from which an observer model 
reconstructs the full state. This paper shows that using sec- 
ond order observers is a viable technique for reconstructing 
the unmeasured states of structures under mildly restrictive 
conditions. Moreover, the computational advantages of the 
second order observer as compared to a first order observer 
indicate that significantly larger observer models may be uti- 
lized. Numerical examples are used to demonstrate the per- 
formance of second order observers. The implications of sec- 
ond order observers in the development of Control/Structures 
Interaction (CSI) technology is discussed. 

1. Introduction 

Recent emphasis on the optimal interdisciplinary design 
of actively controlled spacecraft has led to the formation 
of a focused research activity within NASA called Con- 
trol /Structures Interaction (CSI). The aim of the CSI ac- 
tivity is to coordinate University. Industry and Government 
research to develop the technology required to make active 
control of flexible spacecraft routine. To this end. several 
investigators have begun to exploit the second order form of 
the differential equations which describe structural dynam- 
ics. The controllability and observability of linear second 
order models has been studied by Laub and Arnold 1 . Hughes 
and Skelton 1 and by Bender and Laub. 3 

Park and Belvin 4 ’ 3 demonstrated the computational advan- 
tages of simulating the coupled CSI equations in second order 
form. Unlike classical first order simulation of the CSI equa- 
tions. simulation in second order form enables the symmetry 
and sparsity of the structural equations to be exploited. Sig- 
nificant computational efficiency is gained by treating the CSI 
equations in second order form. Moreover, a natural parti- 
tioning of controls and structures is indicated by maintaining 
the structural equations in second order form. 3 

The linear quadratic regulator (LQR) control synthesis tech- 
nique has excellent robustness properties. Safonov and 
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Athens* have proven, when the full state is available, that 
the LQR technique provides 60* of phase margin and infinite 
gain margin. However, for space structures the full state is 
not available and an observer must be used in conjunction 
with LQR. No guarantees can be made on the robustness 
properties for the observer based LQR controller. 7 Thus the 
observer plays a key role in the study of CSI when full state 
feedback is employed. 

Observer state reconstruction suffers from model errors, 
measurement noise and unmodeled dynamics. Unmodeled 
dynamics result from the use of reduced order models. 3 To 
reduce the level of unmodeled dynamics, the size of the ob- 
server could be increased, provided the observer solution can 
still be performed in real-time. The computational advan- 
tages associated with the second order observer may permit 
larger observers with no penalty in solution time as compared 
to the first order observer. 

Second order observer models have received recent attention 
in the literature. Hashemipour and Laub* utilized an optimal 
observer known as the Kalman Filter for discretized (time and 
space) second order structures models. In addition, robust 
computational procedures for solution of the Kalman Filter 
estimation error covariance matrices have been developed for 
second order models in Ref. 10. Unfortunately, the continu- 
ous time Kalman Filter has only been successfully developed 
in first order form. 

It seems paradoxical that the optimal way to estimate the 
states of a second order system is with a first order observer 
model. Although there is more design freedom when using 
a first order model, the disadvantage of numerically simu- 
lating the CSI equations in first order form may prove that 
sub-optimal second order observers are most practical, ft 
is from this perspective that some analytical and numerical 
results are presented herein to highlight the advantages and 
disadvantages of using second order observers to estimate 
the states of realistic spacecraft. 

In section 2. the governing CSI equations are presented using 
nomenclature similar to Refs. 9 and 10. It includes the clas- 
sical first order observer and a logical form for the second or- 
der observer. Stability and performance of the observers are 
discussed in Section 3. Practical computational aspects of 
the observers are examined in Section 4. Section 5 presents 
numerical results which demonstrate the viability of second 
order observers. In addition, implications of the second order 
observer on CSI technology are discussed. 
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2. Controlled Stru cture Simulation Equations 

A control-structure interaction system can be represented as 
shown in Fig. 1. The linear time-invariant discretized equa- 
tions of motion for the structure (plant) may be described 
by the following set of continuous-time equations: 

Mq + Dq + Kq = f + Bu + Gw 

q(0) = qo. q(0) = q> (2-1) 

z - Hjq + H„q + v 
u = F x q + F a q 

where M, D, K are the mass, damping and stiffness ma- 
trices, respectively. These matrices, typically derived using 
finite element modeling, are symmetric and sparse. The vec- 
tor f represents known, state independent, applied forces. B 
is the actuator location matrix and u is the state dependent 
control force. The matrix G describes the disturbance loca- 
tion and w is the disturbance vector. The displacement and 
velocity initial conditions are given by q, and q,. respectively. 
The measured output of the system is s where v is a vector 
of measurement noise. The matrices Ha and H. represent 
the displacement and velocity sensor locations and F, and 
F 3 are feedback gain matrices of the controller. 

The conventional technique for dealing with multi-variable 
control problems is to convert the system to first order state- 
space form. 

x = Ax + Ef+Bu + Gw (2*2) 

z = Hx + v, x(0) = Xq 
u = Fx 

where 



H = [H„ H.], F = [F l F a ] 



The control gain F may be synthesized by a number of tech- 
niques as given in standard texts such as Ref. 11. When 
full state feedback is used, and the full state cannot be mea- 
sured. an observer model must be developed to estimate the 
states of the entire structure from the measured states . The 
form of the first order and second order observer models is 
discussed below. 

2.1 First Order Observer Model 

The observer model is usually constructed by augmenting 
(2.2) by a state correction term such that 

i = Ax + Ef + Bu + L(z - Hx) (2 • 3) 

x(0) = 0 



Figure 1. Typical Control/Structure Interaction System 



and A. E. B and H are a reduced order description of the 
structural system, actuator and sensor locations. The vector 
x is the estimated state. Adding the state correction term 
L(z — Hx) forces the observer model to track the measured 
states of the structure. Since the disturbance is unknown, it 
is not included in the observer. 

The observer gain matrix L can be synthesized by a num- 
ber of techniques. If covariance matrices of the noise and 
disturbance processes in (2.1) are available, the Kalman Fil- 
ter can be constructed.' -11 The Kalman Filter provides an 
optimal balance of observer performance and noise rejection: 
however, a first order observer is required. Since a structure 
is represented with second order models, it seems natural 
to choose a second order observer model for estimation of 
structural states. The second order observer is described 
below. 

2.2 Second Order Observer Model 

A linear observer model could be derived by adding a state 
correction term to a reduced order form of (2.1). such that 
the estimated states q and q are computed from 

Mq + Dq -I- Kq = Ef + Bu + MLaT 

q(o) — 0, 4(0) =0 (2-4) 

where M. D and K represent a reduced order model of the 
structure and 7 is a state correction term of the form 

7 = z - (H rf q + H.4) 

The second order observer form given by (2.4) may be solved 
efficiently by keeping the control and state dependent state 
correction terms on the right-hand-side (RHS) of the equa- 
tion. This permits the symmetry and sparsity of the observer 
mass, damping and stiffness matrices to be exploited by the 
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computational algorithms. References 4 and 5 present a par- 
titioned solution procedure for this type of equation. 

From (2.4) one observes the absence of Li. which must be 
null for a second order observer to be constructed. The ab- 
sence of Li is needed to establish the identity relation q = q. 
A second order system is written in first order form by us- 
ing the relation that the time derivative of position is equal 
to the veleocity as shown in the top row of equation (2.2). 
To transform a first order system to a second order one. 
the identity relation must exist. The effects of L» = 0 are 
described in the next section. 


3. Stability and Performance of Observers 

To study the ability of the observer to predict the states of 
the structure, the equation for the error between the actual 
states and the predicted states. « = (q-q) is examined. To 
simplify the discussion, we consider an unreduced observer 
A. = A, B as B, H = H. The deterministic error is given 
by 

a f ~ LjHrf I — LjH, 1 f3 . l! 

® ~ l -MK - L a Ha -M D - L?H, J e ' 3 *1 

From (3.1). it is seen that the error asymptotically goes to 
zero, provided the matrix [A - LHj is stable. The next 
section describes the stability limitation of second order ob- 
servers using a simple single-degree-of-freedom model. 


3.1 Observer Stability for a SDOF Model 


A deterministic single-degree-of-freedom (SOOF) model con- 
sisting of a spring (k). a damper (d) and a mass (m) is used 
in this section to highlight differences in stability between the 
first order and second order observer models. The system 
equations in first order form may be written for the structure 
as 

*-[—*%• -J/»]’ I + [l/m] ,+ [i/m! u I 3 ' 2 ) 

and for the observer as 

i = f i — l,h» i . r o i f 

[— fc/m — L?Hj -d/m-L J H.] x+ ll/mJ f 

+ [tym] U + [&£ I*] * ( 3 • 3 ) 

The stability of the observer can be determined from the 
poles of the error equation 


*-[- 


- L1H4 1 — L t H. 

k/m — L3H4 —d/m — LjH, 




(3-4) 


The poles of (A - LH) in (3.4) are 

\ ia = a±jw (3-5) 


where 

— (d/m + L a H, + LiHj) 
a 2 


w 3 = [4fc + 2L |£ ijd 4- 4mL a H^ — rrj(L[Hj + L a H u ) 3 
-2L a H,d — 4L 1 H,fc - d a /m)/4m 

The pole locations given by (3.5) show the effect of the ob- 
server gains on the error between the actual and predicted 
states. The L, gain increases the rate of decay of the error 
when position measurements are used. The L t gain also ef- 
fects the frequency of the error response when velocity mea- 
surements are used. The L* gain uses the velocity mea- 
surement to increase the rate of decay and the L, gain uses 
position measurement to effect the frequency of the error 
response. 

From the above considerations, distinct differences are seen 
in the first order and second order observers (L t = 0). The 
first order observer can produce a stable reconstruction error 
utilizing position and/or velocity measurements. The second 
order observer poles are given by 

-(d/m + LjH.) 

o = ■ - ■' 1 

2 

a 4fc + 4mLjH* — m(LjH,)* - 2L?H.d - <P/m 

u = . 

4m 

which show velocity measurements must be available to aug- 
ment the error decay rate. For example, if the structure was 
undamped, (d = 0). and no velocity measurements were 
made (H.=0) then the first order observer would have the 
real part of the poles located at 

„ -L t Hz 
2 

Thus the first order observer would be stable. However, 
the second order observer would be unstable, a = 0. Even 
though stability is gained by a small amount of damping, 
(d 0). the second order observer would have poor per- 
formance without augmenting the damping through velocity 
feedback. Observer performance is discussed in the next sec- 
tion. 

3,2 Performance Considerations 

The performance of the observer depends on the eigenvalues 
and eigenvectors of the matrix [A-LH]. The eigenvalues de- 
termine the speed and rate of observer convergence. Rule of 
thumb choices for the eigenvalues should be somewhat larger 
than the closed loop system [A + BP], but not exceedingly 
so. to prevent problems due to high frequency noise. The 
eigenvalues and eigenvectors of the observer matrix can be 
placed to maximize noise and disturbance rejection. Because 
the second order observer restricts the design. Lj = 0. some 
of the freedom in placing the eigenvalues and eigenvectors is 
lost. 

Robust methods for designing second order observer gains 
remains a fruitful area of research. The results herein have 
synthesized Kalman Filter gains and subsequently ignored 
the Li gain. This approach may or may not lead to a stable 
observer. If L a is not stablizing it must be augmented. Even 
if L3 is stablizing some loss of performance is expected, par- 
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ticularly in terms of noise rejection. However, this may be 
a reasonable penalty for the increased computational speed 
gained by using a second order observer model which could 
be translated into reducing the level of unmodeled dynam- 
ics. The model error (uncertainty) would be the same for 
the first or second order observer. Future research is needed 
to determine the robustness characteristics of second order 
observers. 


4. Computational R equirements for Simulation 

Solution of the CSI equations is typically performed using 
digital computations. It is this consideration that leads to 
the study of second order observers for CSI technology. As 
shown in the following paragraphs, the second order observer 
provides the potential for larger observers and/or faster com- 
putational solutions than can be realized with the first order 
observer 

4.1 First Order Obser ver Solution 

The first order observer given by (2.3) may be solved by 
several methods. One popular method is to compute the 
matrix exponential of 

A. = (A + BF-LH] (4-1) 

so that 

Sc-* 1 = 4(A)x" + r(A)(Ef + Lx") (4 • 2) 

where 

4(A) m eA’ k , T(A) * / t X -‘ds 

Jo 


and h is the temporal step size. 

Defining N as the number of structural degrees of freedom, 
then 9(h) is a 2 N by 2 N matrix. 

Because 4(A) is computed once for time invariant systems, 
the matrix exponential computation does not significantly im- 
pact the total time required for a long simulation. However, 
computing the matrix exponential for large systems does 
cause accuracy and storage problems. 4(A) occupies 4 N 2 
storage locations since the matrix is typically not symmet- 
ric or sparse. Independent of the number of sensors. (4.2) 
requires, as a minimum, multiplication of 4(A) by a vector 
at each time step. This represents 4N* operations at each 
step of the simulation. 

An alternative approach to the matrix exponential route is to 
discretize the equations of motion in time to derive difference 
equations. This approach has been used by the present au- 
thors in Refs. 4 and 5. Midpoint implicit integration formulas 
were used to discretize (2.3) to derive a linear equation 

gjn+1/a _ gfv+i/a 

x" + ‘ = 2x" +l/a -x" (4-3) 

where 

S = I-Ja„, g = J(Ef" +1/a + Lx" +, / J )+x" 
i 2 


To solve (4.3) an L-U decomposition is performed once for 
time invariant systems and subsequently a lower and upper 
triangular system is solved at each time step. Since the L-U 
decomposition of S is not sparse or symmetric, 4N 7 storage 
locations are required. In addition, the two triangular sys- 
tem solutions require 4 N 2 operations. Thus, the first order 
observer requires the same computational resources for both 
the matrix exponential approach and the time discretization 
approach. 


4.2 Second Order Observer Solution 

The second order observer model given in (2.4) can be solved 
using the time discretization approach as presented in Refs. 
4 and 5. A summary of the procedure is presented herein. 


We express a set of the mid-point implicit formulas with the 
step size A: 


q n+,/a = q" + *4 , 6 = A/2 

<r i/a =<r+«r ,/2 

q n+l = 2q" +l/a - q" 


(4-4) 


The selection of the mid-point implicit formula (or the trape- 
zoidal rule) for controlled structures is due to its minimal 
frequency distortion and no numerical damping characteris- 
tics. 


Implicit time discretization of (2.4) using (4.4) yields the 
following difference equation: 

Sq"* 1 ^ = g" +I/a 
S = M + SI) + S 7 K 

|*+t/a _ + Bu» +1 / a + ML 2 <y n+l / 3 ) 

4-lCf (n n + Sa n \ + «Do" ( 4 ' 5 ) 


q"+» = 2q" +l/a - q" 

4~ ,/a = (*"+!/’-*")/* 


From (4.5) it can be seen that u’ v+1 ^ a and 'y" +1 /* are re- 
quired to numerically solve for q" +1/la . The partitioned so- 
lution procedure 5 ’*, uses equation augmentation to predict 
u "4i/2 an( j y»4i/a thereby maintaining the symmetry and 
sparsity of the S matrix even after L-U decomposition. For 
the second order observer. S is an N by N matrix. If the 
bandwidth of S is proportional to N by a constant a: 

Bandwidth = aN 

then, only aN a storage location are needed. Most impor- 
tantly. only 2 N(aN + 1) operations are needed for the two 
triangular system solutions at each time step. 

Comparing the second order observer minimum solution time 
to the first order observer minimum solution time yields 

(: 2aN * + 2 N) a + j? 

4W* ~ 2 

Note this is an approximate count of the number of oper- 
ations as the right-hand-side of (4.2) and (4.5) require ad- 
ditional operations. Nevertheless, for a typical bandwidth 
of 10 percent, and JV 1, the second order observer re- 
quires approximately 5 percent of the number of operations 
as does the first order observer. (Reduced order observer 
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models based on a modal description of the structure yields 
aN = l!) Such a drastic reduction in computations may 
lead to new applications of observers in CSI technology. 

8. Results and Discussion 

To evaluate the performance of the second order observers 
numerically, several examples are presented. The first exam- 
ple is relatively simple, a beam with an unreduced observer 
model. The second example is more realistic in that a large 
truth model is used in conjunction with a reduced order ob- 
server. The encouraging results from these examples prompt 
us to consider serious use of second order observers in CSI 
technology. 



Properties Open loop frequencies (Hz) 

E - 40E + 6 lb/in. 2 *i - 5 6 - 06 

G - 2.4E + 6 lb/in. 2 h - 233.22 


5.1 Beam Example 


p - 1.1 32E - 4 lb-s 2 /in. 


f 3 - 556.73 


To study the effects of approximating a first order observer 
by a second order observer, the beam model shown in Fig. 2 
has been used for numerical studies. The beam is. subjected 
to a step load for 0.002 seconds and a regulator control sys- 
tem is activated at t=0.01 sec. A finite element model was 
constructed with 8 planer beam elements to form a 27 de- 
gree of freedom representation of the beam. A three mode 
modal space control law was used with the performance con- 
straint that the vibration amplitude must be less than 0.025 
in. within 0.1 sec. after the control system is activated. Lat- 
eral position and velocity measurements were made at each 
of the three locations indicated in Fig. 2 for a total of 6 mea- 
surements. A Kalman Filter was constructed to predict the 
unmeasured states of the system. Results for the unreduced 
structure and observer are shown in Fig. 3. The observer, 
shown by the dashed line, rapidly converges to the true po- 
sition of the beam lateral deflection (measured at the point 
of loading). 


To determine the viability of a second order observer. L t 
was set equal to zero and the observer equation was solved 
in second order form. The response of the second order 
observer is shown in Fig. 4. The second order observer is 
slower to converge to the actual beam deflection than the first 
order observer. Nevertheless, the second order observer did 
converge, thus, the observer stability was maintained when 
L, was neglected. 

The beam response with the control law based on a perfect 
observer (all states are reconstructed with no error), the first 
order observer and the second order observer is presented in 
Fig. 5. Although the second order observer produced no- 
ticeable performance degradation in the initial time after the 
control system was activated, it achieved nearly the same 
level of vibration attenuation after t=0.1 sec. as did the per- 
fect observer. Thus, the second order observer is quite viable 
for this example where both position and velocity measure- 
ments were available. 


f 4 - 2693.6 

Figure 2. Pinned-Pinned Beam Structure 




Time, sec 

Figure 4. Beam and Second Order Observer Response 
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Figure 5. Comparison of Beam Response with Different 
Observers 

S.2 Earth Pointing Satellite Example 

A more realistic example of spacecraft being studied for ac- 
tive vibration control is the Earth Pointing Satellite (EPS) 
shown in Fig. 6. This model is a derivative of a geostationary 
platform designed for study of Earth Observation Sciences 
(EOS). Two flexible antennas are attached to a truss (bus). 
Typical missions involve pointing one antenna to earth, while 
tracking with the other antenna. 

The truss has a three meter square cross section. It con- 
sists of 51 mm diameter by 1.59 mm wall thickness graphite 
epoxy tubes. The modulus of elasticity and mass density of 
the tubes are 275.9 E+09 N/M 9 and 3250.0 Kg/M 3 , respec- 
tively. Although the antennas have the same tube dimensions 
as the truss, nonstructural mass make the effective tube den- 
sity ss 9750.0 Kg/M 3 for the two antennas. The finite ele- 
ment representation of the EPS consisted of 570 degrees of 
freedom. Table 1. lists the first 20 frequencies of the EPS 
structure. 


Table 1. EPS Vibration Frequencies (Hz.) 


1. 

0.000. 

2. 

0.000 

3. 

0.000. 

4. 

0.000 

5. 

0.000. 

6. 

0.000 

7. 

0.242. 

8. 

0.406 

9. 

0.565. 

10. 

0.656 

11. 

0.888. 

12. 

0.888 

13. 

1.438. 

14. 

1.536 

15. 

1.776. 

16. 

1.776 

17. 

3.026. 

18. 

3.026 

19. 

3.513. 

20. 

3.531 



Figure 7. EPS Antenna and Second Order Observer 
Response 


To study the viability of second order observers for this struc- 
ture. sensors were located on the 15 M antenna hoop shown 
in Fig. 6 to measure the position and velocity of the antenna 
hoop. In addition, sensors were located at the center of 
gravity (CG) to measure the spacecraft attitude and attitude 
rate. The number of measurements used in the construcu- 
tion of Kalman Filter gains varied from 10 to 20 to create 
three different reduced order observer models. Again, the 
second order observer was derived by neglecting the gain 
which can result in an unstable observer for some problems. 
The observer was studied independent of a feedback control 
law by simply applying an impulsive loading to the open loop 
structure. An impulse was applied as an initial angular ve- 
locity about the z-axis of 100.0 rad/sec between the 15 M 
antenna and the truss. Note that this excitation will excite 
both rigid body and flexible vibrations o{ the EPS. 

Figure 7 shows the structure and observer displacement re- 
sponse for 15 modes in the observer model, and a full 570 
degree of freedom model of the EPS. Although the observer 
is stable, considerable error exists at the higher frequencies. 
Figure 8 shows the error between the EPS states and the ob- 
server states as the number of modes in the observer model 
was increased from 10 to 20. These results show the im- 
portance of larger observers in state estimation of realistic 
spacecraft. 
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Figure 8. Reconstruction Error Convergence for EPS 
Observers 


0. Concluding Remarks 

Classical multivariable control observer gain synthesis algo- 
rithms produce first order observer models. Results have 
been presented which show that approximating the first or- 
der observers by a second order observer can be successfully 
carried out if velocity measurements are available. 

This paper, in conjunction with previous works by the au- 
thors. has shown that second order observer models are 
computationally advantageous. The increased computational 
speed can be translated to larger observers and/or faster 
sampling rates. Larger observers allow the designer to con- 
trol more ‘modes' of the structure which is increasingly im- 
portant for future spacecraft missions. On the other hand, 
if the designer chooses, he may translate the increased com- 
putational speed into faster measurement/control updates. 
This may prevent phase distortion in some applications and 
may actually enable observer based controllers when large 
bandwidths are necessary. 

Advances in C5I technology will occur by exploiting the phys- 
ical attributes of structures and structural models. This pa- 
per should serve to stimulate interest in the use of second 
order observers for state estimation of structures. Second or- 
der observers may prove sub-optimal as compared to the first 
order observer, nevertheless, practical considerations such as 
computational speed may be the decisive factor in developing 
CSI technology. More research is needed to determine if full 
state feedback is the control law of preference for controlled 
spacecraft. If the answer is affirmative, the second order ob- 
server promises to play an important role. However, much 
research is still needed in the areas of directly synthesizing 
second order observer gains and on determining robustness 
of the second order observer. 
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16. Abstract 

The use of Linear Quadratic Regulator (LQR) control synthesis techniques implies 
the availability of full state feedback. For vibration control of structures, 
usually only a limited number of states are measured from which an observer model 
reconstructs the full state. This paper shows that using second order observers 
is a viable technique for reconstructing the unmeasured states of structures 
under mildly restrictive conditions. Moreover, the computational advantages of 
the second order observer as compared to a first order observer indicate that 
significantly larger observer models may be utilized. Numerical examples are 
used to demonstrate the performance of second order observers. The implications 
of second order observers in the development of Control s-Structures Interaction 
(CSI) technology is discussed. 
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